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Abstract 
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"jrt I In this paper a wave is generated by an initial data whose support is locahzed 

at the outside of unknown obstacles and observed in a limited time on a known 
closed surface or the same position as the support of the initial data. The observed 
data in the latter process are nothing but the back-scattering data. Two types of 
^ ' obstacles are considered. One is obstacles with a dissipative boundary condition 

'^ . which is a generalization of the sound-hard obstacles; another is obstacles with a 

—1. I finite refractive index, so-called, transparent obstacles. For each type of obstacles 

CN ' two formulae which yield explicitly the distance from the support of the initial data 

f->^ . to unknown obstacles are given. 
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1 Introduction 

In [9] the author has introduced a simple method for the reconstruction issue of some 
class of inverse obstacle scattering problems. The observation data used therein are given 
by a wave field observed in limited time on a known closed surface surrounding unknown 
obstacles. The wave is generated by the initial data with compact support outside the 
surface. It is shown that the method yields the distance from a given point outside the 
surface to obstacles provided the obstacles are sound-hard or penetrable. It is a kind of 
an enclosure method since the obstacles can be considered as being contained in an open 
ball centered at oo determined by the distance. 

The enclosure method was introduced in [6, 7] firstly for inverse boundary value prob- 
lems governed by elliptic equations. In [8] the author has found its applications in several 
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inverse initial boundary value problems governed by heat and wave equations in one- 
space dimensional case. Quite recently in [11, 12] the method has been extended to 
three-dimensional problems for the heat equation and in [10] the framework for inverse 
initial boundary value problem governed by the heat equation together with a result when 
the background medium is isotropic and inhomogeneous has been established. 

This paper is a continuation of [9] and the purpose consists of two parts. First we 
extend the range of applications of this new method to more general or another kind of 
obstacles. We consider obstacles with a dissipative boundary condition which is a gener- 
alization of the sound-hard obstacles; obstacles with a finite refractive index, so-called, 
transparent obstacles. Second we consider so-called inverse back-scattering problems for 
those obstacles. Therein the observation data are measured in a limited time at the same 
position as the initial data. Needless to say this is a very important class of inverse 
problems for waves and the solution has many possibilities of applications, for example, 
nondestructive testing such as diffraction tomography, subsurface radar, microwave to- 
mography, ocean acoustics tomography, etc. We show rigorously that the method works 
also for inverse back-scattering problems for obstacles mentioned above. 

1.1 Obstacle with a dissipative boundary condition 

First we consider an inverse obstacle scattering problem which is described by the classical 
wave equation outside an obstacle with a dissipative boundary condition. 

Let D be a nonempty bounded open subset of R^ with smooth boundary such that 
R^ \ D is connected. Let 7 and /3 be functions belonging to L°°{dD) and satisfy 7 > 0. 
Let < T < 00. Given / G L^(R^) with compact support satisfying supp / fl D = let 
u = Uf{x,t) denote the weak solution of the following initial boundary value problem for 
the classical wave equation: 

d'^u - Am = in (R^ \D)x ]0, T[, 

u{x,0) =0 inR^\D, 

dtu{x,0) = f{x) mR^\D, 

dvj 

— - -f(x)dtu - p(x)u = OondDx ]0, T[. 
ou 

Here z/ denotes the outward normal to D on dD. Note that the boundary condition with 
/3 > is called the dissipative boundary condition, however, in this paper we include also 
the case /3 < in this terminology. Note also that a formal computation yields 



:i.i) 



E'(t) = -2 / -f(x)\dtu\^dS <0 

JdD 

E{t)= I _{\dtu\'^ + \Vu\'^)dx+ I /3(x)|Mprf5, tG [0, T]. 
J-r:a\d JdD 

Let VL he a 
R^ \ r2 is connected. In this subsection we consider the following problem: 



IdD 

where 

Eit) = [ _, 

Let Q he a bounded domain of R^ with a smooth boundary such that D G Q and 



Inverse Problem I. Assume that D^ 7 and (5 on dD are unknown. Extract information 
about the location and shape of D from the wave field Uf{x,t) given at all x G dQ and 
t G ]0, T[ for a fixed / G L^(R^) with compact support satisfying supp / fl f2 = 0. 

Thus dQ is a mathematical model of the closed surface on which the wave field is 
measured. 

Let B denote an open ball. In this paper as done in [7] we introduce two conditions 
on initial data /: 

(11) f{x) = a.e. xeR^\B; 

(12) there exists a positive constant C such that f (x) > C a. e. x E B or —f{x) > C 
a.e. X E B. 

Note that (11) and (12) imply that supp/ = B. 
Set 



w{x) = Wf{x, t) = e '^^Uf{x, t)dt, x G R \ D, r > 0. 



;i.2) 



We denote by u also the outward normal to Q on dQ. 

We assume that 

• there exists a positive constant C such that the one of the following two conditions 
is satisfied: 

(Al) 7(x) <l-C' a.e. x G dD; 

(A2) 7(x) > 1 + C" a.e. x G dD. 

The first result of this subsection is the following formula which employs w together 
with dw/du on dfi. 

Theorem 1.1. Let B satisfy B 011 = ^ and f e L'^{Yl?) satisfy both (II) and (12). Let 
V = Vf{- ,t) G i/^(R^) denote the weak solution of the modified Helmholtz equation 



{A-T^)v + f = inR^. 
If the observation time T satisfies 

T > 2dist{D,B) - dist{n,B), 



;i.3) 



;i-4) 



then we have: 

if (Al) is satisfied, then there exists a Tq > such that, for all t > Tq 



T^w - -—V \dS>0; 
an \ou ou 



;i.5) 



if (A2) is satisfied, then then there exists a tq > such that, for all t > tq 



T^w - -—V ]dS <0; 
an \ov ou 



1.6) 



In both cases the formula 



lim — loe 

r^«, 2r 



/ dv dw \ 
an \du du 



-dist{D,B) 



;i.7) 



is valid. 



Remark 1.1. w in Theorem 1.1 is explicitly given by the formula 

1 f e^'^\'-^~y\ 

v{x) = Vf{x,T) = — -. 7f{y)dy (li 

Atx Jb \x — v\ 



and thus computable from / in principle. 

From the data Uf{x,t) on dQx ]0, T[ one can compute w together with dw/du on 
dQ. The computation process makes use of the uniqueness of an initial boundary value 
problem for the wave equation in (R^ \ r2)x]0, T[. We do not repeat this description in 
detail and see [9] for the procedure. Thus formula (1.7) is: a solution to Inverse Problem 
I provided 7 satisfies (Al) or (A2); an extension of Theorem 1.1 in [9] which treated 
sound-hard obstacles, that is, the case when 7 = /3 = on dD. 



Note that dist {D, B) + J\dB\/47r coincides with the distance from the center of B 
to D. Thus by choosing a suitable initial data which is independent of D, from (1.7) one 
gets the distance of -D to a given point in R^ \ Q. 

By Proposition 1.1. in [9] (1.4) ensures that T > l{dB,dD,dQ), where 

l{dB, dD, on) = inf{|x -y\ + \y-z\\xedB,yedD,ze d^}. 

Thus (1.4) does not contradicts the finite propagation property of the signal governed 
by the wave equation. It should be emphasized that we never make use of the finite 
propagation property of the signal governed by the wave equation in any form and this is 
an interesting point of our method since we could find restriction (1.4). 

It should be pointed out also that Theorem 1.1 gives a characterization whether 7 < 1 
or 7 > 1 in terms of a solution of the wave equation over a finite time interval provided 
7 is constant and 7 7^ 1. 

Our method covers also another very important case which employs the data observed 
on the support of the initial data not on dQ, that is, the back-scattering one. More 
precisely we consider the following problem. 

Inverse Problem I'. Assume that D, 7 and f3 on dD are unknown. Let B satisfy 
i? n D = 0. Extract information about the location and shape of D from the wave field 
Uf{x, t) given at all x G -B and t G ]0, T[ for a fixed / G L^(R^) with support B. 

A new finding in the enclosure method what we want to emphasize is the following 
result. 

Theorem 1.2. Let B satisfy BnD = $ and f e L'^(R^) satisfy both (II) and (12). Let 
V = Vf{ ■ ,t) E H^(IV^) be the weak solution of (1-3). IfT satisfies 

T >2dist{D,B), (1.9) 

then we have: 

if (Al) is satisfied, then there exists a tq > such that, for all t > tq 

f{w-v)dx>0; (1.10) 

if (A2) is satisfied, then there exists a tq > such that, for all t > tq 

f{w-v)dx<0; (1.11) 



In both cases the formula 

I r 

-dist{D,B), (1.12) 



lim — lo£ 



f{w — v)dx 

B 



is valid. 



In formula (1.12) we make use of w{x,t) for x E B which can be computed directly 
from Uf{x,t) for {x,t) G Bx ]0, T[ via (1.2). Since supp / = B, this is a back-scattering 
data over finite time interval ]0, T[. Thus (1.12) gives a solution to Inverse Problem I'. 
Note also that B can be arbitrary small. 

Theorem 1.1 ensures that even the case when T is not greater than 2dist {D, B), how- 
ever satisfies (1.4) one can extract dist (-D, B) via formula (1.7). This is an expression 
of an advantage of making use of data Uf on the whole of dQ. In back-scattering case, 
restriction (1.9) on T is quite natural expected or optimal one because of the finite prop- 
agation property of the signal in the wave phenomena, however, note that in the proof 
we never make use of this property in any form. 

In the framework of the Lax-Phiilips scattering theory [13], Majda [14] considered the 
case when /3 = and D is strictly convex. He clarified the leading term of the scattering 
amplitude s{6,u,X) at high frequency which is the Fourier transform of the scattering 
kernel S{s,6,u). However, our data are different from his ones and in this paper we never 
assume that D is strictly convex. 

1.2 Obstacle with a finite refractive index 

Second we consider an inverse scattering problem for an obstacle with a finite refractive 
index. 

Let < T < oo. Given / G L^(R^) with compact support let u = Uf{x,t) denote the 
weak solution of the following initial boundary value problem: 

a{x)d'^u - Au = inR^x]0, T[, 

M(x,0) = OinR^ (1.13) 

dtu{x,0) = f{x) inR^ 

where a is a function belonging to L°°(R^) and satisfies a{x) > C a.e. x G R^ for a 
positive constant C . 
We assume that 

• there exists a nonempty bounded open set D with a smooth boundary such that 
a{x) = 1 a.e. x G R'^ \ -D; 

• there exists a positive constant C such that the one of the following two conditions 
is satisfied: 

(Bl) a{x) <l-C' a.e. x e D; 
(B2) alx) >1 + C' a.e. x e D. 

Thus a has a jump across dD. D is a mathematical model of an obstacle with a finite 
refractive index. 

In this subsection we present a solution to the following problem. 



Inverse Problem II. Let fi be a bounded domain of R^ with a smooth boundary such 
that D G Q and R^ \ f] is connected. Assume that both D and a in D are unknown. 
Extract information about the location and shape of D from the wave field Uf{x,t) given 
at all X G dfl and t g]0, T[ for a fixed / G L^(R^) with compact support satisfying 

supp / n n = 0. 

Set 

rT 

w{x) = Wf{x,T) = / e~^^Uf{x,t)dt, X eR^, T > 0. (1-14) 

•^ U 

Theorem 1.3. Let B satisfy B nU = ^ and f e L'^{Yl?) satisfy both (11) and (12). Let 
V = Vf{ ■ ,t) G H^(R.^) be the weak solution of (1.3). 

IfT satisfy (1.4), then, we have: 

if (Bl) is satisfied, then there exists a tq > such that, for all t > tq 



—w - -—V ]dS>(}; 
an \ov Of 



if (B2) is satisfied, then there exists a tq > such that, for all t >tq 

—w - -—V \dS <0. 
an \ov ov 



In both cases we have 



lim — log 



( dv dw \ ,^ 
an \ov ov 



-dist{D,B). 



;i.l5) 



;i.l6) 



;i.i7) 



The procedure for the computation of both w and dw/dv on dVt from Uf{x,t) given 
at all X G dQ and t G ]0, T[ is the same as Theorem 1.1. Thus this is a solution to Inverse 
Problem II. 

The method works also for the following problem. 
Inverse Problem II'. Assume that both D and a in D are unknown. Let B satisfy 
B r\ D = ^. Extract information about the location and shape of D from the wave field 
Uf{x, t) given at all a; G -B and t G ]0, T[ for a fixed / G L^(R^) with support B. 

The following result gives a solution to Inverse Problem 11'. 

Theorem 1.4. Let B satisfy BnD = ^ and f e L'^{Y{?) satisfy both (II) and (12). Let 
V = Vf{ ■ ,t) E iJ^(R^) be the weak solution of (1.3). 

IfT satisfies (1.9), then we have: 

if (Bl) is satisfied, then there exists a tq > such that, for all t > tq 



f{w — v)dx > 0; 

if (B2) is satisfied, then there exists a tq > such that, for all r > tq 

/ f(w — v)dx < 0. 
Jb 

In both cases we have 



lim — lo£ 

T — 5>oo 2r 



B 



f{w — v)dx 



-dist(D,B). 



Xll 



;i.i9) 



;i-20) 



In the framework of the Lax-Philhps scattering theory there is a resuh by Maj da- 
Taylor [15] for the case when a is smooth on D and a{x) ^ 1. However, the data are 
given by the back-scattering kernel S'(s, — a;,^), w G S*^ and the obtained information is 
the value of the support function of unknown obstacles and thus different from ours. For 
the study of the leading term of the scattering amplitude at high frequency for a strictly 
convex obstacle with a finite refractive index see Majda- Taylor [15] {a{x) > 1 on D) and 
Petkov [16, 17] {a{x) < 1 on D). 

1.3 Construction of the paper 

A brief outline of this paper is as follows. Theorems 1.1 and 1.2 are proved in Section 
2. In Subsection 2 first we formulate what we mean by the weak solution of (1.1). It is 
based on the notion of the weak solution in [1] and an application of the existence theory 
therein. We see that w given by (1.2) satisfies the modified Helmholtz equation with an 
unknown inhomogeneous term outside D and boundary data on dD in a weak sense. In 
Subsection 2.2 two expressions for the integral 

(%-§^.)rfS (L21) 

are established. One is called the local expression which yields a bound of the absolute 
value of (1.21). Another is called the global expression and yields the lower bound of 
(1.21) in the case when (Al) is satisfied: upper bound of that in the case when (A2) is 
satisfied. Those bounds are presented in Subsections 2.4 and 2.5, respectively. Theorem 

1.1 is a direct consequence of those bounds. The proof of Theorem 1.2 is described in 
Subsection 3.6. It is based on the following asymptotic formula which connects two data 
in Theorems 1.1 and 1.2 for an arbitrary fixed T < oo: 

_f(w~v)dx= f (^w-^v]dS + 0(T~^e-^^). (1.22) 

R3\D Jan \ou av ) 

Using this together with the obtained bounds in Subsections 2.4 and 2.5, one gets imme- 
diately the conclusion of Theorem 1.2. 

Section 3 is devoted to the proof of Theorems 1.3 and 1.4. The order of the presentation 
is parallel to that of Section 2. Starting with the formulation of the weak solution of 
(1.13), we present two expressions of (1.21) for w given by (1.14). Both expressions are 
of global type since they involve integrals over the whole space. Using those expressions, 
we give a bound of the absolute value of(1.21) in Subsections 4.2; the upper/lower bound 
of (1.21) for case (B1)/(B2). Theorem 1.3 is a direct consequence of those bounds and 
an asymptotic formula which connects two data in Theorems 1.3 and 1.4 similar to (1.22) 
enables us to obtain all the conclusions of Theorem 1.4. 

Section 4 consists of some remarks and open problems. In Subsection 4.1 we show 
that the case 7 = 1 is exceptional in one-space dimensional case. Moreover, the complete 
asymptotic expansion of (1.21) in one space-dimensional case is presented. In Subsection 

4.2 corresponding to results in one-space dimensional case we propose some open problems 
and point out related results in some references together with future direction of our 
method. 



2 Proof of Theorems 1.1 and 1.2 

2.1 The weak solution and the governing equation for Wf 

The contents of this subsection almost parallel to the corresponding parts of Section 2 in 
[9] which is the case when 7 = /3 = on dD. 

We write u' instead of dfU. It has been observed in [9] that for the enclosure method 
the notion of the weak solution described on pp 552-66 in [1] is appropriate since we 
consider only the solution on a finite time interval. 

Set V = H\R^ \ D) and H = L\R^ \ D). Applying Theorem 1 on p 558 in [1] to 
(1.1), we know that: given u^ E V and u^ E H there exists a unique u satisfying 

u e L\0, T; V), u' e L2(0, T; V), u" G L2(0, T; V) 

such that, for all G K 

< u"{t), (f) > +a{u{t), 0) + b{u'{t), 0) = 0, a.e. t g]0, T[, (2.1) 

and m(0) = m" and u'{0) = u^, where 

a{u,v) = / _\/u -Vvdx + / jSuvdS, u,vgV 

JR3\D JdD 

and 

b(u, v) = 'juvdS, u,v E V. 

JdD 

Note that this is the case when / = in their notation and from their proof one knows 
that the restriction (5.8) on p 553 for b is redundant in this case, that means that 60 
therein can be identically zero and b = bi. Note that at this step we do not make use of 
condition 7 > on dD. 

In this section we say that this u = Uf for u^ = and -u^ = / is the weak solution of 
(1.1). Then w = Wf given by (1.2) belongs to V; it follows from integration by parts and 
(2.1) multiplied by e""^* that w satisfies for all G V 

/ _\/w -Vddx + / c(x,T)w(f)dS + / _{t'^w — f)(hdx 

J-R?\D JdD Jn?\D 

(2.2) 

-tT 



f _F(/)dx-e-^^ f G(l)dS, 

JR3\_D JdD 



— e _ 

where 

c{x,t) = 7(a:)r + l3{x), 

F{x, t) = u'{x, T) + Tu{x, T), 

G{x) = 7(x)-u(x, T). 
This means that, in a weak sense, w satisfies 

(A - r2)«; + / = e-^^F{x, r) in R^ \ D, 
— — = dx, t)w + e~'^'^G(x) on dD. 

Of 



In particular, from (2.2) for G C^(R^ \ D), we have (A — t'^)w + f{x) = e '^'^F{x,t) 
in R^ \ D in the sense of distribution and hence Aw G L^(R^ \ D). This yields w G 
Hf^^{R^ \ D) and (A - t^)w + f{x) = e'^^F^x, r) a.e. a; G R^ \ D. Thus we can define 
dw/du\Q^ as VwloQ ■ u G H^^'^{dVL), where Vw\sii is the trace of Vw onto dVl ([4]). 

2.2 Local and global expressions of (1.21) 

In the following proposition we do not assume that supp f Cl D = (/} nor supp / fl ^2 = 0. 
Proposition 2.1. We have two expressions: 

( dv dw \ .^ 
an you ou J 

,- c{x,T)v\wdS — I (w — v)fdx (2-3) 

dD \du J Jn\W ' 

""^ / GvdS-e-^^ I _Fvdx; 

JdD Jn\D 



f dv dw \ ,^ 
an \ou ou 



f (\Vv\^ + T^\v\^)dx- f c{x,T)\v\^dS-e~^^ f GvdS 

Jd JdD JdD 



_\\V(w-v)\^ + T^\w-v\'^]dx+ c{x,T)\w-v\'^dS (2.4) 

R3\D '^ J JdD 

e-^^ f _ F{w -v)dx- e~"^ / _ Fvdx + e""^ / G{w - v)dS 

Jr^kd Jn\D JdD 



(w — v)fdx — / fvdx. 
n\D Jd 

Proof. First we give a proof of (2.3). Let (/? G H^{Q \ D) satisfy (/? = on dVt in the sense 
of the trace. Since the zero extension to R^ \ D of this ip belongs to V ([4]), we have from 
(2.2) 

0= / _Vw -V^pdx + / c{x,T)w(pdS 
Jn\D JdD 

(2.5) 
+ / _{r^w - / + e-^^F)ipdx + e""^ / GipdS. 

Jn\D JdD 

Choose X ^ C'^(R^) such that x = 1 in a neighbourhood of dQ and x{^) = in a 
neighbourhood of D. Since (1 — x)'^\n\D vanishes in a neighbourhood of dfl, (2.5) for 
<^ = (1 - x)v\n\D is valid. ^ ^ 

On the other hand, since xw vanishes in a neighbourhood of D and w G H? (R^ \ D), 



integration by parts yields 

dw 



I _Vw ■V{x'v)dx = / ——xvdS— / _(Aw)xvdx 
Jn\D Jan ou Jn\D 

dw 



f ^vdS - f _{r^w -f + e-^^F)xvdx, 
Jan ou Jn\D 



that is 

dw 



f -^vdS = ( _Ww ■ V(xv)dx + f (t'^w - / + e'^'^F)xvdx. 
Jan ou Jn\D Jq\d 



in\D Jn\D 

From this and (2.5) for (/? = (1 — x)'^\n\D "which satisfies ip = v on dD we obtain 

dw 

'jdb = _ 
ln\D 



/ ——vdS= / _\/w -Vvdx + / c(x,T)wvdS 
Jan ou Jq\d Jao 

+ / S'T-'^w - / + e-^^F)vdx + / e-^^GvdS. 
Ja.\D Jan 



(2.6) 



By the trace theorem one can choose w G H^{Q) such that w = w in Q\D. Note that 
w = w on dQ and dD in the sense of the trace. Since v G H^{Q), we have 

dv f dv 

-—wdS = / -—wdS 
an du Jan du 



= / Avwdx + / Vw ■ Vwdx 
Jn Jn 

- / _(r^f — f)wdx + / _\/v ■ Vwdx 
Jn\D Jn\D 

+ / (r^f — f)wdx + / Vf ■ Vwdx. 
Jd Jd 



On the other hand we have 

dv 



/ Vf ■ Vwdx = / —-wdS — \ (Av)wdx 
Jd JdD du Jd 

= / —-wdS — / (r^f — f)wdx, 
JaD du Jd 



that is, 

/ —-wdS = / (r^t> — f)wdx + / Vf ■ Vwrfx. 
JaD du Jd Jd 

Therefore, we obtain 

f dv f r f dv 

/ ^—wdS= / _(T^v-f)wdx+ / _Vv ■Vwdx+ / ^—wdS. (2.7) 

J9Q az/ Jn\D Jn\D JdD du 

A combination of (2.6) and (2.7) gives (2.3). 
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Next we give a proof of (2.4). Write 

/ —wdS= —{w-v)dS+ —vdS. 2.8 

JdD dv JdD OU JdD OU 

By the trace theorem given G V^ one can choose G i/^(R^) such that = 0|r3\-d- 
Since f G H'^{D) and Af — r^f + / = a.e. x E D and = on 9-D in the sense of the 
trace, integration by parts yields 

/ -^(t)dS = f {Vv- Vcpdx + Th(p)dx - I f(pdx. (2.9) 

JdD OU Jd Jd 

On the other hand v satisfies 

- / Vt; ■ V^dx -r^ f .^ v<^dx = - I ^ f^dx, Vy? G /J^(R^). 

Substituting ^p = 4> into this identity and dividing R^ = D U (R^ \D), we obtain 



_(Vv ■ V0 + T^v(f))dx = / _ f(l)dx - / (Vv ■ V0 + r^vcfydx + / fcjydx. 
R3\D jr3\d jd Jd 



A combination of this and (2.9) gives, for all G V^ 

dv 

(pdti + / 

Combining this with (2.2), we obtain 



/ -^(t)dS+l _(Vt;-V0 + r\0)rfx= / _/0rfx. (2.10) 

JdD OU JR3\D Jn^XD 



f l^-c(x,T)v](l)dS=f c(x,T)e(l)dS + f _ | Ve ■ V0 + r^e0| rfx 

JdD \0U J JdD JR3\D '^ J 

+e-"^ / _F(f)dx + e-^^ I G(t)dS, 

where e = w — v. This means that e satisfies, in a weak sense 

(A - r2)e = e-^^F in R^ \ D, 

— cix, r)e = — 7" — h c(a;, t)v + e~'^'^G on 9D. 

(7Z/ (7Z/ 



(2.11) 



Substituting e for in (2.11), we obtain 

dv 

edt) = I c(a;, r)(|e|^ + t>e)d^' + / 

i\D 



f ^edS=f c{x,r){\e\^ + ve)dS+ f {\Ve\^ + r^\e\^)dx 

JdD du JdD JRJiXD 

+e-^^ f _Fedx + e-^^ f GedS 

JR3\D JdD 



(2.12) 



Now from this together with (2,3), (2.8) and the identities ve = —\v\'^ + vw and 

dv 



f ^vdS= f {\Vvf + T^\vf)dx- I fvdx, (2.13) 

JdD du Jd Jd 
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we obtain (2.4). 
D 

Now assume that supp/ fl fi = 0. From (2.3) and (2.4) we have the following two 
expressions which we call the local expression and global expression, respectively: 



f dv dw \ ,^ 
an \ou ou 



/ \ — cix, t)v wdS — e '^^ GvdS — e '^^ _ Fvdx; 

J an you J Jao Jn\D 



(2.14) 



f dv dw \ ,^ 
an \ az/ ou 



(2.15) 



f {\Vv\^ + T^\v\^)dx- f c{x,T)\v\^dS-e-^^ f GvdS 
JD JdD Jan 

\\\/(w — v)\'^ + t'^Iw — vl"^} dx + / c(x,t)\w — vl'^dS 

R8\D 1^ J JdD 

+e-^^ f _ F{w - v)dx - e~"^ / _ Fvdx + e""^ / G{w - v)dS. 
Jk^\d Jn\D Jao 

For convenience we set 

We make use of (2.14) to give an estimation of |/(t)| from above and (2.15) of /(r) from 
below/above when (A1)/(A2) is satisfied. This is the role of (2.14) and (2.15). 

2.3 An estimate of |/(r)| from above 

In this subsection we derive the following estimate as r — > oo: 

g2rdist{D,B)|j(^)|^ 0(^3/2)^ (2.16) 

Lemma 2.1. Let e = w — v. As r — > oo we have 

l|e|lL^(R3\D) = Oie-'^^ + e~'^dzst(D,B)^ (2.17) 

and 

l|Ve|li.(R3\i5) = 0{r\e-'^^ + e-'^d^stiD,B)^y (2.18) 

Proof. Rewrite (1.12) as 

/ _i\^e\^ + T^\e\^)dx = - f c(x,T)(\e\^ + ve)dS+ f ^edS 

J-R^XD JdD JaD OU 



-e-^^ 



[ _Fedx+ f Geds] 

^JR3\D JdD J 

12 



(2.19) 



We have 



t Ov 

c(x,T)(\e\'^ + ve)dS + / 7— erfS* 

dD JdD OU 



-T / 7Me + ve 

IdD 



)dS+ I {^-(3\e\^-(5vt\dS. 

JdD \OV J 



Completing the square and 7 > on dD give 



■/ 7(|e|2 + w)rf5=- / 7e + !J dS + j 
JdD JdD 2 Ja 



IdD 4' 



\v\ dS 



< I ^Ivl'^dS. 

dD 4 



Moreover we have 



/ (^e-/3\e\'-/3ve]dS<C f (\e\^ + \v\'+^']dS, 

JdD you J JdD \ OU I 



-tT 



GedS 



and 



-tT 



Fedx 



dD 



1 



<- / |e|2rf5 + e-2-^/ \G\^dS 

- 2 \JdD ' ' i9D ' ' 



From (2.19)-(2.24) we obtain 



- 2 I ' Jn^\D ' ' ' 



(2.20) 



(2.21) 



(2.22) 



(2.23) 



/ iFpdx I , Vr7 > 0. (2.24) 

Jl{?\D J 



\Ve\^dx+ r^-^] I \e\^dx 



R3\D 



2 / jR-i\D 



<{C + l]f \e\^dS+f (c + ^)\v\^dS + C 



dD 



2 J JdD 

2 XJdD 2t]^ Jr3\d y 



dD 



dv 



dv 



dS (2.25) 



Here we cite the well known inequality (e.g., [4]) that there exists a positive constant 
K = KiVt \ D) such that, for all z e H^{n\D) 

f \z\^dS <Klr]^ f _\Vz\^dx + r]-^ f _\z\^dx] , \/r]>0. (2.26) 

JdD \ Jn\D Jn\D J 
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Applying (2.26) to the first term in the right-hand side of (2.25), we have 

n2" 



|Vepdx+ r2-!L) / \e\^dx 



R»\Q 



2 / jR-i\n 



{l - K,' (C . 1) } l_ IV.P.. ^{(r'-'j)- Krr' (c + i) ) j^^_ \.H. 



dD 



dv 



du 



-2tT 



dS + ^^if \G\^dS+^f \F\^dx 
2 \JdD 2if Jr3\d ) 



(2.27) 
Now choosing a small r^, from (1.8) and (2.27) we obtain (2.17) and (2.18) as r — > oo. 

D 



Now from (2.26) for ri = 1/^/t, (2.17) and (2.18) we obtain 



l^l|2 nf^f^.-'^T'T I ^-2rdist (D,B)\\ 



This together with 11^11^2(9^) = 0{e-^^^^^^^'^'>) yields 



I'"^IIl2(9D) 



0(r(e-2-^ + e^2rdist(D,B))) 



(2.28) 
(2.29) 



Since 



dv 



dv 



O^ g-rdist(D,B). 



L2{dD) 



it follows from (2.29) that 



-—wdS 

dD OV 



+ 



dD 



c{x,T)vwdS 



Moreover we have 

and since 
we obtain 



GvdS 



(-)(^^3/2(^g-2rdist(D,B) ^ ^-^T{T+(\\si {D,B)-^y (^2. 30) 
Q(g-r(T+dist(D,B))) ^2.31) 



dD 



U,|| _ _ /nfp-Tdist(f7,B)N 



-,-rT 



Fvdx 



n\D 



(-,/^g-r(T+dist(Q,B))N 



(2.32) 



Now from (2.14), (2.30), (2.31) and (2.32) we obtain 

|J(^)| ^ (-,(^^3/2(^g-2rdist(D,B) _^ g--r(T+dist (D, B))^) _^ ^g-r(T+dist (Q, B))^) 

and thus 

g2rdist(D,B)|J^^^^|| ^ Qj^^3/2j^]^ _^g-r(T-dist(Z),i?))^| ^ ^g-r(T-2dist (D, i?)+dist (C, B))^|_ (^2. 33) 

Write 

T - dist {D, B) = (T - 2dist (D, B) + dist (f], 5)) + (dist {D, B) - dist {Vl, B)). 

This together with dist (£*, 5) > dist (^2, B) and (1.4) gives T > dist (D, B). Now from 
this together with (1.4) and (2.33) we obtain (2.16). 

Remark 2.1. From this last part we know that, for the proof of (2.16) it suffice to assume 
only T > 2dist (-D, B) — dist (f], B) instead of stricter condition (1.4). 
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2.4 Case (Al). An estimate of /(r) from below 

In this subsection we prove that: there exists a positive constant fi such that 

hm inf r^e'"^^^* (^' ^^/(r) > (2.34) 

T >-00 

provided 7 satisfies (Al). 

Having (2.16) and (2.34), we immediately obtain (1.5) and (1.7). 

For the proof of (2.34) first we study the asymptotic behaviour of the following integral 

J{t)= f {\Vv\^ + T^\v\^)dx - f c{x,T)\v\^dS^. 
Jd JdD 

Since B Cl D = (l), f{x) = a.e. x & D and thus (2.13) gives another expression: 

J{t) = I {^-c{x,T)v\vdS. (2.35) 

Lemma 2.2. There exists a positive number ^ such that 

lim r^e^rdist(D,B)j^. > q_ ^2.36) 

T S-OO 

Proof. Since 

ou An Jb |x — i/p An Jb \x — y\'^ 

we have 

'^(^) = TT^ dS^ k{x, y) - -l{x, y) jrf{y)f{y')dydy' 

[Ati]^ JdD JbxbK r } \x — y\\x — y'\ 

where 

, , , (y — x) ■ v(x) , , 

\x-y\ 

Ha;,!/) = — . f^ P(a:). 

F ~ ^1 

We divide the integrand of Jir) into two parts. Set d = dist {dD, B) and Ai = 
{{x, y) G dD X B\\x — y\ = d}. It is easy to see that d = dist (D, B). 

In what follows we denote by Br{z) the open ball centered at a point z with radius 
R. Given 5 > define 

W5 = U(,„,j,,)e>,(aD n 55(^0)) X (fin 5^(1/0)) X (5 n ^^(yo)). 

The set W5 is open in dD x B x B and contains the set of all (x, y, y) with (x, y) G Al. 
Here we state the following two claims concerning the Ws- Their proofs are almost same 
as those of Claims 1 and 2 in [12]. Just rewrite the proofs by replacing D and dQ in the 
previous definitions of Ai, Ws and F{x,y,y') in [12] with dD and B; {x — y')/\x — y'\ 
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in the previous definition of F[x,y,y') in [12] witli z/(x). By tliis reason we omit tlie 

description. 

Claim 1. Given e > tfiere exists a ^i > sucli tliat for all (x, y, y') G Wsi it holds that 

(y — x) ■ uix) , , , , ,, , 

^ — - — ^ > 1 - e, x-y <d + e, x-y' <d + e. 
\x-y\ 

Claim 2. Given 5i > there exists a ^2 > such that if {x, y, y') E dD x B x B \ Ws^, 
then \x — y\ + \x — y'\ >2d + 62- 

Let C be the constant in (Al). Give e = C'/2 in Claim 1 and choose Si in Claim 1. 
Next choose 62 in Claim 2 corresponding to 61 already chosen. 

By claim 1, we have, for all {x,y,y') G Ws^ k{x,y) > 1 - e - (1 - C") = C'/2, 
\x — y\ < d + e and \x — y'\ < d + e; by claim 2, we have e"''*^''^"^'^'^"^''^ < e^'^'^'^e''^'^^'^ 
if {x,y,y') G dD x B x B\ Wsi- These together with estimate f{y)f{y') > C^ for a.e. 
y E B gives 

'^(^) = TT^ / dS^dydy' k{x, y) - -l{x, y) jrf{y)f{y') 

[4:TrY Jws^ V r / \x — y\\x — y'\ 

+ 7r^ / dS^dydy' k{x, y) - -/(x, y) ^ jrf{y)f{y') (2-37) 

(47r)"' J9DxBxB\Wi^ V T J \x - y\\x - y'\ 

{C'/2 - T-^Ci)tC^ f ^^ M.^..'^'r{\x-y\ + \x-y'\) , ^^^^-2rd^-T52\ 



> ^ ' w.' r / dS.dydy'e'-^y^-y^^^-'y'^^ + 0(. 



where 

^' = dist(A5) "^ l'^ll^°°(^^)- 

Choose (xo,yo) ^ -^- Since \xo — yo\ = d, we have 

\x — y\ + \x — y'\ < 2\x — Xo\ + \xq — y\ + \xq — y'\ 

< 2\x - xo\ + 2d+ \yo - y\ + \yo - y'\. 
From the definition of W^^ we have 

dSJydy'e'^^^''~y^+^'"y'^^ > f dS^ f dy f rf^'g-d^-J/l+l^-y'l) 

Ws^ JdDnBs-^{xo) JBnBs-^iyo) J BnBg^iyo) 



,5i(a;o) yjBfiBs^iyo) 



IdDnBs^ (xo) \JBnBs^ (yo) J 

(2.38) 

Here we make use of two facts which are essentially same as Claims 3 and 4 in [12]: 
for all S > we have 



iBnBsiyo) 

and 



liminfr^ / e""l^°~^'rfy > 

liminfr^ / e-=^"l^-^»lrf5^ > 0. 

T — ^oo JdDnBxixn) 



IdDnBsixo) 
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From these, (2.37) and (2.38) we see that (2.36) for /x = 7 is true. 



n 



It follows from (2.26) that 

/ _\Ve\'^dx + T^ f _|ePrfx+ / /3\e\'^dS 

J-R»\D JR3\Z) JdD 



> (1 - K\\^\\L^(^aD)V) / _ iVeI'dx + / _ \Ve\'dx 



+ {r' - K\\/3h^^eD)V'') f _ M'dx + r' f _ 

Jn\D JR3\n 



> (r' - K\\(3\\L^^9n)V-n -M'dx 

JR3\D 



> — / lel'^dx 
2 Jr3\d 



el'^dx 



(2.39) 



for a fixed r] with 1 > i^||/3||L°°(9D)^ and r > J2K\\l3\\L°^(Qr))r] . Moreover we have 



, _\e\^dx + e ^'^ I _Fedx = -, _ 
2 7r3\d J-r?\d 2 Jr3\d 



= -rT 



re 



2 e"^'^'^ 

dx — 



2t^ 7r3\£) 



iFrdx 



-2tT 



> 



\Frdx 



and 



2t2 Jr3\d 
r / 7|e|2rf^ + e-"^ / GedS 

JdD JdD 



(2.40) 



r / 7 

IdD 



-tT 



e + ^^u{x,T) 



2 



rf5 



4r JaD 



7|M(x,T)prfa; (2.41) 



-2tT 



> 



'y\u{x,T)\'^dx. 



At JdD 

Note that we have made use of the assumption 7 > on dD. Now it follows from (2.15), 
(2.39), (2.40) and (2.41) that 

nr) = J(r)+f \Ve?dx + T^ I lel^dx + f c(x,T)\e\^dS 
Jr3\d Jr3\d JdD 



+e 



-tT 



Fedx — e 



-tT 



R3\D 



n\D 



Fvdx + e-^' I GedS 

dD 



-2tT 



-2tT 



>J{r)-%r^l \F\^dx-—^ — / -f\u(x,T)\^dx-e~^^ f Fvdx 
2r^ Jr3\z) 4r JdD Jn\D 

= J{t) + 0(e-2-^ + re-^(^+dist (n,B))y 
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and thus 



,2rdist {D, 



^2rdist (D, B) 



^t^^ZTLilfil(lJ,B)j/\ y ^fj.^2Tmiil(V,B)j/\ 



-2T{T-dist{D,B)) I fi+1 ~T{T~2dist {D, B)+d[st (n, B))-s 



Now from this together with (2.36) we obtain (2.34). 



2.5 Case (A2). An estimate of /(r) from above 

In this subsection we prove that there exists a positive constant /x such that 

hminff-r''e2-^i^*(^'^)/(r)'l >0 



provided 7 satisfies (A2). 

Having (2.16) and (2.42), we immediately obtain (1.6) and (1.7). 
The proof of (2.42) starts with rewriting (2.12) as 



dD 4c(x, 



r 



' dv 
dv 



cix, r)v \ — e ^ G 



2 -2tT 

dS 



At"^ yR3\D 



\F\^dx 



(2.42) 



C(X,T 



dD 



e-:.^l(?^-c(x,r).]-e--G' 



2c{x,t) Wdiy 



dS 



(2.43) 



Since 



and 



iVeP + r^ 



R3\D 



-tT 



2r2 



dx. 



-tT 



'-^^^ 



1, ,0 e 



-2tT 



At* 



' dv 



2c(x, r) \\dv 



c{x,t)v\ -e~^'G 



>^|e| - 



2' ' 4c(x,r)2 
it follows from (2.43) that 



— c{x, t)v — e '^^G 



dD 2 7r3\d \ 2 / 



< 



dD 2c{x, t) 



'dv 
dv 



dx, t)v — e "^ G 



2 -2tT 

dS + 



2r2 Jn.3\D 



iFl^dx 
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and thus 



c(x,T)\e\^dS+ / (\Ve\^ + T^\e\^] dx 

dD Jr^\D ^ 



< 



dD C[X,T 



'|-c(x,rJ-e--G 



2 .-2rT 

dS + ^ 



\F\^dx 



dD C[X,T) 



dv 
du 



— ex, TjV 






It follows from (2.17) that 

e-^^ / _Fedx = 0(r(e-2-^ + e-(^+dist(D,B)))) 

and from (2.28) that 

e-^^ I GedS = 0(ri/2(g-2.T ^ ^-.(r+dist (d,b)))^ 

Now applying (2.31), (2.32), (2.44), (2.45) and (2.46) to (2.15) we obtain 

Hr) < J{t) + 0(r(e-2-^ + e--(^+dist (d,b)))) ^ Q(^g-r(T+dist (n.B))^ 
where 



J{r) = J{t) + 



9D C{X, T) 



dv 
dv 



c(x, t)v 



dS. 



From (2.35) one gets 



J{r) 



dD C{X,T) \Oh' J OU 



T 



v 



1 



-dSr 



X / ( k{x, y') l{x, y') ) (k{x, y) /(x, y) 



in J JdD c{x, r 

X_ \ p-T-ilx-y'l+l^-yl) 



(2.44) 



(2.45) 



^2.46) 



(2.47) 



\x — y \\x — y\ 



fiy')fiy)dy'dy, 



where k is the same as that of Lemma 2.2 and 

(y-x)- iy{x) 



k{x,y) 



l{x, y) 



\x-y\ 
(y-x) ■ iy{x) 



\x-y\ 



From (A2) we have k{x,y') < —C and thus 



k{x,y') - -l{x,y') < -^, {x,y') e dD x B 

T I 
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for all T > tq with tq » 1. Now give e = 1/2 in Claim I of Lemma 2.2 and choose 
6i and corresponding 62 in Claim 2. It is easy to see that one can apply the same 
argument as done in Lemma 2.2 to J{t) where k{x,y) — {1/T)l{x,y) plays the same role 
as k{x,y) — (l/r)/(a:,t/) in Lemma 2.2. Thus we conclude that there exists /i > such 
that 

liminf|-rV^dis*(^'^)j(r)l>0 

and thus from (2.47) one gets (2.42). 

2.6 Proof of Theorem 1.2 

Since i? fl D = 0, one can find a bounded open set Q with a smooth boundary such that 
R^ \ r2 is connected, D G il and B (IQ = ^. We consider I{t) for this Q. 
Since supp / fl f] = 0, it follows from (2.3) that 



dv 

wdS — I c{x,T)wvdS — e '^ / GvdS — e "- _ 

ln\D 

Substituting (p = w into (2.10), we have 

dv 

wdS = I 

and thus (2.48) becomes 



I(t)= f T^wdS- f c(x,T)wvdS-e-^'^ f GvdS - e'^^ f _Fvdx. (2.48) 

JdD OU JdD JdD Jn\D 



/ —-wdS= / _fwdx— / _CVv -Vw + T'^vw)dx 

JdD OU Jr3\D JR3\D 



Ht~) = / f(w — v)dx—< / (Vw -Vv + (t'^w — f)v)dx + / c(a;, r)wf rfS* > 

Jr.«\D [Jll«\D JdD J 

-e-^^ f GvdS-e-^^ f _Fvdx. 

JdD Jn\D 

Now it follows from this and (2.2) for = f |J^3^-J) that 

/(r) = / _f{w-v)dx + e-^^ f _Fvdx. (2.49) 

Jn^XD jR3\n 

Since ||F||i2(R3\;D) = 0(r) and 

||^||L2(R3)<r-l/|U2(R3), (2.50) 

(2.49) yields 

/ _f(w-v)dx = I(T) + 0(T-^e-^^). (2.51) 

JR3\D 

Now (1.10) and (1.12) follow from (2.51) together with (2.16) and (2.34) in case (Al); 
(1.11) and (1.12) follow from (2.51) together with (2.16) and (2.42) in case (A2). 



20 



3 Proof of Theorems 1.3 and 1.4 

Set V = H\R^) and H = L'^{R^). Applying Theorem 1 on p 558 in [1] to (1.13), we 
know that: given u^ & V and u^ E H there exists a unique u satisfying 

MeL2(0, T-V), u'eL\0, T;V), d/dt{C{u'{-)) e L\0, T;V') 

such that, for all G K 

< —C{u'{t)), <p > +a{u{t), 0) = 0, a.e. t e]0, T[, (3.1) 

(JjL 

and u{0) = vP and 'u'(O) = -u^, where 

aiu.v) = I Vu -Vvdx, u,v &V 

and C : H — > H is the bounded linear operator defined by 

C{u) = au, u E H. 

Note that this C satisfies (5.11) on p 553 in [1] under the condition a{x) > C a.e. x G R^ 
for some positive constant. Note that at this step we do not make use of condition 
a{x) = 1 a.e. x G R^ \ -D. 

Note also that this is the case when 6 = in Theorem 1. However, the equation is 
homogeneous, that is, / = in their notation and by virtue of this their proof covers also 
this case. 

In this section we say that this u for u^ = and -U"^ = / is the weak solution of (1.13). 

We see that w given by (1.14) belongs to V and applying integration by parts to (3.1) 
multiplied by e~'^'^, we obtain, for all G V^ 

/ (Vw ■ V0 + T^aw(f))dx - I af(t)dx + / F(t)dx = 0, (3.2) 

■/R3 JR? JR3 

where 

F{x) = e-^'^a{x){u'{x,T) + Tu{x,T)). 

This means that w is the weak solution of the following equation 

(A - aT^)w + af = F in R^ (3.3) 

By the same reason as w in Section 2 we have w G H? (R^) and (3.3) holds a.e. x G R^. 
Since v satisfies (1.3) in the weak sense, we see that w — v satisfies, in a weak sense, 

(A - ar^) {w-v) = {a- 1) {t% - f) + F mR\ (3.4) 

3.1 Two expressions 

We start with the following two key expressions. Note that in the proposition we do not 
assume that supp / fl f2 = 0. 
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Proposition 3.1. We have: 



( dv dw \ .^ 

dCl \0U OP 



= I {|V(w -v)\^ + r^a|w - t;p} dx + t'^ f {I - a)\v\'^dx 

- / (1 - a)nw - v)dx +f Fiw- v)dx + f fiav - w)dx - / Fvdx; 
JR3 JR3 Jn Jn 



(3.5) 



(3.6) 



f dw dv \ ,„ 
—-V - -—w dS 
an you ou J 

= I \\W(v-w)? + T'^\v-w?]dx + T'^ I {a-l)\w?dx 

+ / {a — l)f{v — w)dx — / F{v — w)dx + / f{w — av)dx + / Fvdx. 
JR3 JR3 Jn Jn 

Proof. The proof of (3.5) and (3.6) are rather simpler than that of (2.3) and (2.4). First 
we give a proof of (3.5). Since v G H^{Q) and satisfies (A — r^)f + / = a.e. x E Q, 
integration by parts gives 

/ —-wdS = / (r^f — f)wdx + / Vw ■ Vwdx 
Jan ov Jn Jn 

and by the same reason for w we have 

/ -T—vdS = / Vw ■ Vvdx + t"^ awvdx — / a fvdx + / Fvdx. 
Jan dv Jn Jn Jn Jn 

These yield 

/ T— w — T^v \ dS = t'^ (I — a)wvdx + / fiav — w)dx — / Fvdx. (3.7) 

Jan \ov ov ) Jn Jn Jn 

Write the integral in the first term of the right hand-side as 

{1 — a)wvdx = / {1 — a)v{w — v)dx + / {l — a)\v\^dx. (3.8) 

n Jn Jn 

On the other hand, from the weak form of (3.4) and integration by parts we have 

V iw — v)\^ dx — t"^ I a\w — v\^dx 
R3 7r3 



that is, 



{a — \)v(w — v)dx — I {a — l)f{w — v)dx + / F{w — v)dx, 

R3 JR3 Jr3 



r / (1 — a)v(vu — v)dx 

JR3 



/ {\V(w -v)]"^ + T'^a\w -v\'^}dx + (1 - a)f(w -v)dx + F(w-v)dx. 

JR3 Jr3 Jr3 

(3.9) 
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Since 1 — a{x) = outside D, from (3.7), (3.8) and (3.9) we obtain the desired expression. 
Second we give a proof of (3.6). 
We change the role of v and w. Introducing 

f = f--F, F = --F, 
a a 

one can rewrite equations (1.4) and (3.3) as 

{A-T^)v + f = F inR^ 



and 

Then v — w satisfies 



{A - aT^)w + af = inR^ 

(A -T^){v-w) = F-f-{A- T^)w 

= F-f-{A-aT^ + {a- l^jw 

= F - f + af + {1 - a)T^w 

= {l-a){r^w-f) + F. 



Integration by parts gives 

dw 



and 



/ ——vdS = / (Aw)vdx + / Vw ■ Vvdx 
Jan ov Jn Jn 

= / ar'^wvdx — / afvdx + / \/w ■ Vvdx 
Jn Jn Jn 

= / a{T^w — f)vdx + / Vw ■ Vvdx; 
Jn Jn 

/ -—wdS = / {Av)wdx + / Vf ■ Vwdx 
Jan ou Jn Jn 

= / {t'^v - / + F)wdx + Vv ■ Vwdx 
Jn Jn 

I Vf ■ Vwdx + t"^ I vwdx — / fwdx + / Fwdx. 
Jn Jn Jn Jn 



From these we obtain 



( dw dv \ ^^ 
—-V - -—w dS 
an \ov ou 



T^ {a — l)wvdx + / f{w — av)dx — \ Fwdx. 
Jn Jn Jn 



(3.10) 



(3.11) 
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Here rewrite the integral of the first term in the right-hand side: 

{a — l)wvdx = / (a — l)w{v — w)dx + {a — l)\w\'^dx. (3-12) 

It follows from the weak form of (3.10) that 

\\/(v — w)\'^dx — t"^ / \v — wl'^dx 

R3 Jr3 

(1 — a){T'^w — f){v — w)dx + / F{v — w)dx. 

R3 iR3 



This yields 



r' I (a — l)w(v — w)dx 

/r3 



{\\/{v -w)\^ + T^\v -w\^}dx- {1 - a)f{v -w)dx+ F{v-w)dx 

R3 JR» JR3 



{\V(v - w)r + r> - wr}dx + / {F+(a- l)f}(v - w)dx. 
R3 jr» 

Now from (3.11), (3.12) and (3.13) we obtain 

f dw dv \ ,„ 
-—V - —-W dS 
an \ov ov 



(3.13) 



Since 



= T^ {a- l)\wfdx + / {\V{v -w)\^ + T^\v - wf}dx 

+ / {F + {a — l)f}{v — w)dx + / f{w — av)dx— / Fwdx. 
JR3 Jn Jn 

F + {a-l)f = -F + {a-l)f, 

this is nothing but (3.6). 
D 

Now assume that supp / fl f] = 0. Set 

Since a{x) = 1 a. e. a; G R^ \ D, (3.5) and (3.6) become the following two expressions of 

Kir): 

K{t) = I ^ {|V(w -v)\^ + T^a\w - t;p} dx 

(3.14) 
+r^ / (1 - a)|t;p(ix+ / F{w -v)dx- / Ft;(ia;; 

JD JR» Jf7 
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-K( 



T 



— / F{v — w)dx + / Fvdx. 
JR3 Jn 



(3.15) 



3.2 An estimate of |A'(r)| from above 

We describe the proof of the following estimate as r — > cxd: 

g2.dist(D,B)|^(^)|^Q(^2)^ 

Since {a{x) — l)f{x) = for a.e. x G R^ it follows from (3.9) that 

/ I |V(w — f )p + r^Q;|w — f l^j- (ix 



r / (1 — Q;)f (w — v)dx — / i^(t(7 — v)dx. 

JD JR» 



(3.16) 



(3.17) 



Applying a much simpler argument than that for the derivation of (2.17) and (2.18) in 
Lemma 2.1 to (3.17) and using the estimate 



P||2 _ ni''r2p-2^^^ 



we obtain, as r — > oo 



\w - t;||i2(R3) = 0(e-2-^ + e-2-dist (i^.s)^ 



(3.1^ 



|V(w-t;)||i2(R3) 



and thus one gets 



R3 



-F(w — v)dx 



0(r2(e-2^^ + e-2-di^*(^'^))) 



0(r(e-2-^ + e-^(^+dis^(^'^)))). 



(3.19) 



(3.20) 



Since Ht'lUacn) = 0{e'^^^^^^^'^^), from (3.18) one gets also 



Fvdx 



Q(^^^-r(T+dist{n,B)), 



(3.21) 



Q^^-rdist{D,B)-^ and (3.18) to (3.21) that, as 



Now it follows from (3.14), ||f Hizp) 
T — y oo 

\K{t)\ = 0{T\e-^^^ + e-2-dist(D,B))) ^ ^g-r(T+dist (n,B))> 

From this we obtain (3.16). 
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3.3 Case (Bl). An estimate of K{t) from below 

In this subsection we consider the case when a{x) < 1 — C a.e. x & D and prove that: 

hminfrV^dis*(^'^);^(r)>0. 



(3.22) 



Using the completing square 

r^Q;|w — f p + F{w — v) = r'^a 
from (3.14), (3.18) and (3.21) we obtain 



{w — v) + 



2r2 



a 



4r2 



a 



K{r) > C'r^ I |t;|2cix + 0(e-2-^ + re-^(^+dist(^'^»). 

JD 

Here we cite the following fact established in [9]: 

hminfr^e^-di^t^^'^) / \v\^dx>Q. 
This together with (3.23) yields (3.22). 



(3.23) 



(3.24) 



3.4 Case (B2). An estimate of K{t) from above 

Next consider the case when a{x) > 1 + C" a.e. x E D. In this subsection we prove the 
following estimate: 

liminf ('-r^e2-dist(i),B);^(^)\ > q. (3.25) 

We make use of the following completing the square: 

I lO / \iiD iiD iiD 

\v — w\ + [a — l)\w\ = \v\ — 2vw + a\w\ 



aw 



q; — 1 



v\\ 



a 



This together with (3.15), (3.20) and (3.21) gives 

a — 1 , 



> 



C'r 



'^2 



|«||l°°(r3) 



-K{t) >t^ I \v\'^dx - I F{v - w)dx + / Fvdx 

JD a Jn^ Jn 



I |t;prfa; + 0{r{e-''^^ + e'-^^+^ist (d,b)))) ^ ^^^^-.(T+dist in,B))y 

JD 



Now a combination of this and (3.24) yields (3.25). 

It is easy to see that (1.15)/(1.16) and (1.17) follow from the combination of (3.16) 
and (3.22)/(3.25) in the case when (B1)/(B2) is satisfied. 
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3.5 Proof of Theorem 1.4. 

Since B (1 D = (/}, one can find a bounded open set Q with a smooth boundary such that 
R^ \ n is connected, D cQ and 5 n H = 0. 

We make use of (3.7). Since supp / fl f] = 0, (3.7) becomes 

/ I T^w - -TT^A dS = T^ / (1 - a)wvdx - / Fvdx. (3.26) 

Jan you op j Jq Jn 

One can rewrite the first term of this right-hand side as 



r / (1 — a)wvdx = r / (1 — a)wvdx 
Jn JR3 



= r^ / vwdx — T^ awvdx 
Jr3 Jr3 

= — / Vf ■ Vwfia; + / fwdx — t^ awvdx 
7r3 Jr3 Jr3 

= — / (Vw ■ Vf + T'^awv)dx + / fwrfx 
Jr3 Jr3 

= — ( / a fvdx — / Fvdx j + / fwdx 
= / f (i(7 — af)(ia: + / Fvdx 

JR3 Jr3 

/ f(w — v)dx + / (1 — a)/t;(ia; + / Fvdx. 

JR3 JR3 Jr3 



Since (1 — a)/ = in R^, from this and (3.26) we obtain 

/ T— w — -— f rfS* = / f(w — v)dx + / _Fvdx. 
Jan you ov ) jr3 jR3\n 

This together with (2.50) and (3.18) gives 

/ f{w-v)dx=l (^w-^v\dS + 0{r-^e-^^). 
Jr3 Jqq ydu ov ) 

Now (1.18)/(1.19) and (1.20) follows from this together with (3.16) and (3.22)/(3.25) 
in the case (B1)/(B2). 

4 Remarks and open problems 

4.1 One space dimensional case 

Let us consider the meaning of conditions (Al) and (A2) in one space dimensional case. 
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Let < T < oo and a > 0. Given / G -^^(R) with supp / C ] — oo, a[ let m = Uf{x, t) 
be the weak solution of the following initial boundary value problem: 

d'^u - (9^M = in] - oo, a[x ]0, T[, 
u{x,Q) = 0, dtu{x,0) = f{x) in] — oo, a[, (4.1) 

-d^^u{a, t) - -fdtu{a, t) - [5u{a, t) = 0, < t < T, 

where 7 and /3 G R are constant and 7 > 0. This is the case when D =]a, 00 [. We take 
Vt =]0, oo[ and thus dVt = {0}. 

Hereafter we choose initial data / in such a way that supp / C ] — 00, 0[. 

Define 

rT 

w = Wf{x,T) = e '^*Uf{x, t)dt, a; G ] — 00, a[, r > 0. 

J (j 

w" — t'^w + f{x) = e~'^'^F(x, r) in] — 00, a[ 

w'{a) + c{T)w{a) + e'^'^g^r) = 0, 

F{x, r) = dtu{x, T) + tu{x^ T), 

g{T) =-fu{a,T), 

c{t) = '^t + (5. 
Assume that the support of / is compact and satisfies 

supp / C [6 — e, 6] = i? 

with 6 < and e > 0. 

Note that / G i^^(R) with supp / c] — 00, h] always satisfies 



w satisfies 
and 
where 



Now define the indicator function by 

/(r) = -v'{0, t)w{0, t) + w'{0, t)v{0, r), r > 0, 
where 

1 foo _ I _ 

V = v{x, r) = — / e ^1^ y^f{y)dy. 

ZT J-oo 

V belongs to iJ^(R) and satisfies the equation v" — r^f + f{x) = in R. 
We have the following asymptotic formula /(r) as r — > 00. 

Lemma 4.1. Let f satisfy (4-5). As r — > 00 we have 



J(^) = 1 {l-'^)r + /3 ^_2rdist(D,B) 
^ ' 2r(7 + l)r + /3 



-T{h~y) 



f{y)dy 



(4.2) 
(4.3) 

(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4i 



+0(e"^^'^+^^^^^^'-^)^ + 7-e-^(^+^^5^(f^'^))^ 
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Proof. Using one space-dimensional version of (2.14), (2.31) and (2.32), we have 

/(r) = -{v\a) + c{T)v{a))w{a) + 0(e-^(^+di«*(^'^)) + re--(T+dist(n,B))y (4.9) 

Note that dist {D,B) = a-b and dist (fi, B) = -b. 
Since 

this together with (4.5) and (4.7) gives 

v{a) = ^ f e-^'^^-y^f{y)dy, 
2r Jb-e 

v'{a) = -I f e-<'^-y^f{y)dy 

Z Jb—e 

and thus 

v'{a) + c{T)v{a) = i (^ - 1) e-<'^-'^ f^ ^e'^^'-y^ f{y)dy. (4.10) 

Next we compute w{a). Since w G L^(— oo,a) one can write 

1 r 



W{X, T) 



r e-^\^-y\{f{y)-e~^^F{y,T))dy + Ae^^^~^\ (4.11) 

Zt J-oo 



where A is an undetermined constant. From this we have 



w'{x,t) = -- 



(£^ e--(--^)(/(y) - e-^^F{y, r))dy - £ e'^^y-^\f{y) - e'^'^Fiy, r))dy^ 



Substituting (4.11) and (4.12) into (4.3), one gets 

1 /•" 



(4.12) 



e-<''-y\f{y)-e-^''F{y,r))dy + rA 

Z J—QO 

Mr) (^ £^ e~^^''-y\fiy) - e-^^Fiy, T))dy + a) + e-^^gir) = 



From this we obtain 



A = -^^4 V R^ - 1 ) / " e--^^'y^fiy)dy 

2(c(T)+r) \ T ' ' ^y; y 



+ ^^i;L^f£M_lU"e-.<.-».f (,.,),, ^-'^ 



2(c(r) + t) \ T J J-oo c{t) + r 

and (4.11) yields 

w{a,T) 



9{r) 



CiT] + T J-oo 



<''-y^f{y)dy - ^^^ ( r e--("-^)F(i/, r)rfi/ + g{r] 

C[T) + r \J-oo 



(4.13) 
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From (4.4), we see that the integral of the second term of (4.13) has the bound 0{t^^'^). 
From this together with the trivial equality 



e-^i'^-y)f(y)dy = e-^"-'"'^ / e'^^^'-y'^ f{y)dy 

oo J —oo 



we obtain 

w{a, t) = ^^ r e-^^'-y^f{y)dy + 0{e-^^T-^/^). (4.14) 

C{t) + T J-oo 

Now a combination of (4.5), (4.6), (4.9), (4.10) and (4.14) gives (4.8). 
D 

Here we introduce another important assumption on / which means, implicitly, that 
a; = 6 is not a zero point of f{x) with infinite order: 



3/x G R liminf r 



r ¥oo 



b 



b-e 



e-^^'-y^f\y)dy 



> 0. (4.15) 



Having (4.8), one can easily obtain the following theorem. 

Theorem 4.1. Let f satisfy (4.5) and (4.15). Let T > 2dist{D, B) - dist{n, B). 
(i) Assume that 7 7^ 1. We have: 
if '-f < 1, then there exists a tq > such that, for all t > tq 

/(r) > 0; 
if '-f > 1, then there exists a tq > such that, for all t > tq 

I{t) < 0. 
In both cases we have 

lim — log|/(r)| = -dist{D, B). (4.16) 

(a) Assume that 7=1 and /3 7^ 0. We have: 

if /3 < 0, then there exists a Tq > such that, for all t > tq 

/(r) > 0; 

if /3 > 0, then there exists a tq > such that, for all r > tq 

I{t) < 0. 

In both cases we have (4-16). 
Similarly to (2.51), one has 



B 



f{w-v)dy = I{T) + 0{T''e'^'). (4.17) 



This connects two types of data asymptotically. 

A combination of (4.8) and (4.17) yields the following theorem for back-scattering 
case. 
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Theorem 4.2. Let f satisfy (4.5) and (4.15). Let T > 2dist{D, B). 
(i) Assume that 7 7^ 1. We have: 
if '-f < 1, then there exists a Tq > such that, for all r > Tq 

f{w - v)dy > 0; 

B 

if '-f > 1, then there exists a tq > such that, for all t > Tq 

f{w — v)dy < 0. 

JB 

In both cases we have 



B 



1 
lim — loe 



f(w - v)dy 



-dtst{D,B). (4.18) 



(a) Assume that 7 = 1 and /3 7^ 0. We have: 

if P < 0, then there exists a Tq > such that, for all t > Tq 



f{w - v)dy > 0; 
if /3 > 0, then there exists a tq > such that, for all r > tq 

f{w — v)dy < 0. 



B 



In both cases we have (4. 18). 

Thus it is interesting to consider what happens in the case when 7 = 1 and /3 = 0. 

In this case from (4.8) we see that e^"^ *^^''^^/(r) is exponentially decaying as r — )■ 
00 provided T > 2dist (D, B) — dist (f], B). Thus it seems that one has to study more 
about the remainder term hidden in (4.8) by 0(e"^(^+'^^*^*(^'^)) + re"^(^+^^*^*(^'-^))). 
However, in fact, we see that u{x,t) and thus /(r) are independent of a by the following 
argument. 

Given / G L'^(R) with compact support define 

u{x,t) = l r'f{y)dy. (4.19) 

Z Jx~t 

We see that 

dfU + dxU = f{x + t). 

Thus if supp / c] — 00, a], then /(a + t) = for any t > and thus u satisfies —d^u^a, t) — 
dtu{a,t) = 0, t > 0. Of course u satisfies the wave equation in the whole space and the 
initial conditions u{x,0) = 0, ^^^(a;, 0) = f{x). Therefore by the uniqueness of the weak 
solution of (4.1) we have u{x,t) = Uf{x,t) for (x, t) G ] — 00, a[x]0, T[. However it is easy 
to see that the right-hand side of (4.19) is independent of a provided supp / c] — 00, 0[. 

Summing up, we conclude: 

• if 7 = 1 and /3 = 0, then D is invisible by the data Uf{0, t) for < t < T or Uf{x, t) 
for (x, t) E Bx ]0, T[ for any T and / with sup / C ] — 00, 0[. 

Finally we note that one can also determine 7 and /3 after having known dist {D, B). 
More precisely we have the following formula. 
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Theorem 4.3. Let f satisfy (4-5) and (4-15). 

(i) IfT > 2dist{D, B) — dist{Q, B), then as r — > oo we have the following complete 
asymptotic expansion: 



^2Tdist(D,B) 



Kr) 



7-1 1 



/3 



'T{b-y) 



f{y)dy 



2(7+l)r (7 + 1) 



1)2 S 



/3 



n=0 



7 + 1 / r 



n+2 



(4.20) 



(a) IfT> 2dist{D, B), then as r — > oo the function 

^2rdzstiD,B) f f^^_^^^y 



B 



^-rib-y) 



f{y)dy 



has the same asymptotic expansion as (4-20). 

Proof. (4.15) means that there exist a positive constant C and tq > such that for all 

T > To 

(4.21) 



g-T(fe-J/) 

Jb-e 

This together with (4.8) gives 

g2rdist(D,B)j/ N 



f{y)dy > Cr-^. 



-T(b~y) 



f{y)dy 



1 (7-l)r + /3 
2r(7 + l)r + /3 



,Q/2ii-T(T-dist{D,B)) I 2^4+1 -T(T-2dist(i:), B)+dist (f7,B))\ 

Thus expanding the first term of this right-hand side as r — > oo, we obtain (4.20). 

(ii) is a direct consequence of (4.17), (4.21) and (i). 
D 

Note that from the coefficients of the first and second terms of the right-hand side in 
(4.20) one gets 7 and successively /3, 

4.2 Some open problems 

• To study the asymptotic behaviour of (1.21) as r — )■ 00 as done in Theorem 4.3 in 
one-space dimensional case. We think that in three dimensions some restriction of the 
geometry about D should be imposed and have to use (2.14) or (2.52) instead of (2.15). 
For the purpose the asymptotic profile of w on dD or B has to be clarified as done in one 
space dimensional case, see (4.14) 

• The case when 7 = 1 and /3 = seems pathological, however, from a mathematical point 
of view it would be interesting to study the asymptotic behaviour of I{t) as r — )■ 00 in 
three-dimensions since in one-space dimensional case one can not extract the distance of 
S to D. 
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When the data is given by the back-scattering kernel in the Lax-Philhps scatter- 
ing theory, a corresponding positive result in three dimensions has been announced in 
Georgiev-Arnaoudov [2] . 

Note also that Hansen [5] considered a corresponding problem in the Lax-Phillips 
scattering theory for a transparent obstacle with a smooth a on D whose trace onto dD 
is 1, however, its normal derivative on dD is nonzero everywhere. This is the case when 
both of (Bl) and (B2) are not satisfied. He showed that the back-scattering kernel still 
catches the values of the support function of D as the leading singularity. What can one 
say in our problem setting? 

• Can one say something when the data are measured on another fixed ball B' with 
B r\D = ^ and B 0-8 = 0? For example, what happens on the asymptotic behaviour of 
the following integral involving / and another function /' with supp f = B as r — > oo: 



B' 



f'{wf-Vf)dx. 



We think that this is a model of the case when the emitter and receivers of the signal are 
placed on different positions and what we can use is just a single pair of the emitter and 
receiver. 

• It would be interesting to test the performance of formulae (1.7), (1.12), (1.17) an (1.20) 
by using numerically simulated data. 

• Expand the range of applications of the method presented in [9] and here to other 
prototype inverse obstacle scattering problems for acoustic/electromagentic/elastic waves 
or their couplings, etc. In the framework of the Lax-Phillips scattring theory there is a 
result by Georgiev [3] for a moving obstacle. How can one treate the case in our framework? 
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